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Abstract. Time-dependent linear coupling between macroscopic quantum resonator 
modes generates both a parametric amplification also known as a "squeezing 
operation" and a beam splitter operation, analogous to quantum optical systems. 
These operations, when applied properly, can robustly generate entanglement and 
squeezing for the quantum resonator modes. Here, we present such coupling schemes 
between a nanomechanical resonator and a superconducting electrical resonator using 
applied microwave voltages as well as between two superconducting lumped-element 
electrical resonators using a r.f. SQUID-mediated tunable coupler. By calculating the 
logarithmic negativity of the partially transposed density matrix, we quantitatively 
study the entanglement generated at finite temperatures. We also show that 
characterization of the nanomechanical resonator state after the quantum operations 
can be achieved by detecting the electrical resonator only. Thus, one of the electrical 
resonator modes can act as a probe to measure the entanglement of the coupled systems 
and the degree of squeezing for the other resonator mode. 
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1. Introduction 

The quantum behavior of macroscopic resonators has only recently been demonstrated 
in solid-state systems and is currently under intensive exploration. Electrical resonators 
such as a superconducting transmission line resonator coupled to a single Cooper pair 
box or phase qubits and a d.c. SQUID resonator interacting with a flux qubit have shown 
features of single microwave quanta [TJ [21 [3j H] . Nanomechanical resonators capacitively 
coupled to a single electron transistor have been measured, approaching the quantum 
limit with femtometer displacements [5], El [TJ [8] . Recently, a transmission line resonator 
coupled to a nanomechanical resonator has been used to cool the nanomechanical motion 
to hundreds of quanta pj. Besides being a wonderful testing ground for exploring 
quantum physics at the macroscopic level, these systems can provide high-Q harmonic- 
oscillator networks for quantum engineering and quantum information processing |10j . 

Micro-fabricated resonators interact by electromagnetic forces. By controlling the 
circuit parameters, various time-dependent couplings can be generated to manipulate 
the quantum state of the coupled system. In a previous work [IT] , one of us (L.T.) 
studied the effective amplification of the coupling amplitude between a nanomechanical 
resonator and a qubit by parametrically pumping the qubit with a fast pulse. This 
scheme can be applied to produce entanglement between resonator modes as well as to 
generate Schrodinger cat states in the nanomechanical system. In a recent work [12] , it 
was shown that by controlling the coupling between a nanomechanical resonator and an 
ancilla qubit, it is possible to engineer an arbitrary Hamiltonian for the nanomechanical 
resonator mode. In another work [13], a pulse technique is used to form arbitrary 
quantum states of a nanomechanical resonator in the coupled resonator-qubit system. 
In addition, experimental realization of parametric coupling has been achieved for two 
coupled flux qubits [T3] . 

In this paper, we show that quantum features such as entanglement and squeezing 
can be generated and tested using the Gaussian states of macroscopic resonator 
modes coupled in a tunable way. Parametric modulation of the coupling strength 
results in a squeezing operation or a beam splitter operation, necessary ingredients for 
engineering quantum features. Starting from a separable initial state, inseparability (or 
entanglement) can be generated by linear operations. Two physical systems are studied: 
in one system a nanomechanical resonator is capacitively coupled to a superconducting 
electrical resonator mode, in another system two lumped-element electrical resonators 
are coupled to each other through a tunable mutual inductance. In the first system, if 
we start with both the nanomechanical and electrical resonator modes in their ground 
states, careful control over squeezing and beam splitter operations can produce squeezed 
states in the nanomechanical mode. Note that previous work on creating squeezed 
states [15] in nanomechanical-resonator systems uses quantum reservoir engineering, 
feedback control techniques, and an (effective) nonlinear coupling [161 CEO EH fl9] . 
One advantage of the scheme described here is that linear coupling between the solid- 
state resonators, realizable with current technology, can be easily controlled in both 
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magnitude and frequency. In addition, our approach not only provides a method for 
generating entanglement and squeezed states, but it also provides a practical way for 
detecting the squeezing. We show that for Gaussian states, complete information about 
the coupled resonators can be obtained by only measuring the quadrature variances 
of the electrical resonator mode [22] which couples more strongly to the detector than 
the mechanical mode. In this way, we can avoid any difficulties associated with directly 
measuring the nanomechanical mode with high resolution. Here, the electrical resonator 
plays the role of a knob that controls the behavior of the nanomechanical resonator while 
also acting as a detector to probe this behavior. Using a r.f. SQUID-mediated tunable 
coupler, this work can be extended in a straightforward way to a system with two 
coupled electrical resonator modes. For typical dilution refrigerator temperatures (~ 30 
niK), we can pick both resonator frequencies (~ 10 GHz) to operate above the quantum 
limit where the two electrical resonators will definitely be cooled to their ground state, 
making it easy to prepare these modes in squeezed states by linear operations. Also, 
the two modes can be measured simultaneously, which provides a direct observation of 
the cross correlations between the two modes for a clear comparison with theoretical 
predictions. 

Entanglement is the key component for many quantum information protocols such 
as quantum teleportation. Previous work [23] studied the quantum teleportation of 
nanomechanical modes in a purely solid-state network, where high fidelity could be 
achieved for the final state with the assistance of a highly entangled two-mode squeezed 
vacuum state even at finite resonator temperatures. The schemes studied in this paper 
provide a detailed account of how to engineer and evaluate this entanglement to ensure 
the success of the quantum teleportation protocol, as well as to generate one important 
quantum state - a squeezed state of the resonators. 

2. Coupled macroscopic resonators 

In this section, we investigate the realization of linear, tunable coupling between 
macroscopic resonator modes and the parametric modulation of such coupling using 
a circuit architecture. 

2.1. a nanomechanical mode coupled to transmission line resonator 

First, consider a nanomechanical resonator capacitively coupled to a resonant 
superconducting transmission line. Quantum behavior of high-Q electrical modes has 
been demonstrated using superconducting transmission lines [5j . The eigenmodes of a 
one dimensional transmission line resonator of length L between (— L/2, L/2) can be 
obtained by considering the charge distribution along the transmission line, 
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where q(x',t) is the linear charge density at the location x. The lowest even mode has 
the voltage distribution 
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after the quantization of the variables. Here, the frequency of this mode is ufy = 
27r/Lv^c, with I the inductance per unit length and c the capacitance per unit length 
of the transmission line. 

The nanomechanical resonator is located at the middle of the superconducting 
transmission line near x = which is a voltage antinote of this even mode (figured]). 
The nanomechanical resonator is coupled capacitively via a displacement dependent 
capacitance C x . To lowest order, C x = C°(l + x a /d ) depends linearly on the 
displacement coordinate of the nanomechanical motion. The coupled interaction can 
be derived as 

Hint = \cl{l + j)(V x -V{Q)Y (2) 

where the nanomechanical resonator is biased at V x and do is an effective distance 
between the two resonator electrodes. We study the behavior of the lowest 
nanomechanical mode with a displacement x = 5xo(d + a'), where u a is the frequency 
of the mechanical mode, m is the effective mass and 5xq = ^h/2muj a . The interaction 
can be derived as H int = —■^C®V x V(0), assuming that the size of the nanomechanical 
resonator is much smaller than the length of the transmission line. If we apply a voltage 
V x (t) = 2V X sinudt with a drive frequency Ud, the linear interaction can be written as 
Hint = — 2A smudt(a + cfi)(b + b>) with a coupling strength 



(3) 

By adjusting the driving frequency of the voltage, it is possible to generate various 
linear operations, as discussed in more detail in section [31 Hence, the modulation of 
the coupling strength provides an effective tool for controlling the entanglement of the 
coupled resonators. 



2.2. two coupled lumped- element resonator modes 

Alternatively, two electrical resonators can be coupled together using a r.f. SQUID- 
mediated tunable coupler. The circuit schematic (figure [TJ shows a r.f. SQUID placed 
between two lumped-element LC-resonators. The r.f. SQUID acts like an inductive 
transformer where the Josephson junction provides a tunable inductance Lj(5) = 
$ /(27r/ o cos(<5)) and 5 is the phase difference across the junction, I Q is the junction 
critical current, and $o is a flux quantum. The use of small Josephson junctions with 
high current density allows us to ignore the self capacitance Cj of the junction so that 
we remain in an operation regime where uj 2 LjCj 1 for all the relevant frequencies to. 
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Figure 1. (a) A nanomechanical resonator interacts via a capacitance C x with a 
transmission line resonator, (b) Two lumped-clement electrical resonators interact via 
a r.f SQUID-mediated tunable coupler. 



Furthermore, by careful design and operation, we can avoid any direct coupling between 
the two resonators so that a single effective mutual inductance, 

M 2 



MefM 



(4) 



L c + Lj{8) 

derives only from the mutual inductive coupling M between each resonator and the 
central r.f. SQUID transformer coil with geometrical self inductance L c . This leads to 
a coupled interaction of strength [23], 



H int (6) = -M eff (S) 



fljjj a HuJi, 



(a + a+)(S + St), 



(5) 
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where u a and Uf, are frequencies of the two resonators respectively. 

In order to modulate the coupling between the resonators, we must modulate the 
phase difference 8. This can be achieved using a r.f. bias coil which applies an external 
flux cf) x = 27r$ x ./$ a to the r.f. SQUID coil, thereby inducing a current through the 
Josephson junction and thus changing the phase difference 8. The phase difference 8 is 
related to the external flux <p x through flux quantization such that (3 sin(5) + 8 — <p x = 0, 
where (3 = L c /Lj(0). This places some constraints on the design and subsequent 
operation of the tunable coupler [21]. Namely, we would like to operate in a regime where 
(3 < 1 so that the relationship between 8 and cj) x remains non-hysteretic and the inductive 
coupling passes from anti- ferromagnetic (AF) through zero to ferromagnetic (FM). If 
we choose a particular d.c. flux offset and a relatively small amplitude for the r.f. flux 



6 




□ 34 0.36 38 0.4 42 0.2 0.4 0.6 0.8 1 



A/2* f ( au -> 

Figure 2. (a) M e ff/M versus external <P x /2tt, where the (red online) highlight is the 
region where the coupling is modulated. The inset shows the dependence of S on <f) x 
for (3 = 0.85 (b) M e ///M (points) versus time. The thin (red online) line is a sine 
function. 



modulation, we can achieve a roughly linear relationship between the drive flux <f> x and 
the effective mutual inductance M e // without any residual direct coupling. An example 
is shown (figure [2]) for (3 = 0.85 and other device parameters specified later. Thus, 
for a sinusoidal flux drive <f> x (t) = 4>dc — <firf sm(uJdt) at frequency uj, we can produce 
a sufficient parametric modulation of the coupling strength (jSJ) and hence perform the 
squeezing or beam splitter operations. With (pdc = 2.4178 and <p r f = 0.0540, we can 
closely approximate a sinusoidal variation in the size of the effective mutual inductance 
(figure [2]). Slight deviations from the ideal behavior will only appear as negligibly small 
residual direct coupling or small amplitude components at higher harmonic frequencies 
which do not corrupt the squeezing or beam splitter operations. 

A coupling magnitude similar to (jHJ) can be derived for this circuit with 

AM e ff I hw a hw h 

Ao = ^V 2X2^ (6) 

where AM e // is the amplitude of the modulating effective mutual inductance. 

Coupling electrical resonators has several advantages: 1) the coupling between 
the electrical resonators can be made stronger fairly easily through simple design 
modifications, more so than for coupling between a mechanical resonator and an 
electrical resonator, 2) the requirements on device temperature for clearly operating 
in the quantum regime have also been achieved jU [2j [3], 3) electrical resonators have 
already demonstrated sufficiently (high Q's) long coherence times [H El [3], making is 
possible to study their dynamic behavior, and 4) direct measurement of both resonators 
provides more information on the overall quantum behavior of the coupled system. 
Thus, even if pursuing the nanomechanical approach proves to be difficult, it will be very 
promising to pursue coupled electrical resonators in order to investigate the generation 
of entanglement and other quantum physics. 
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3. Linear operations by parametric coupling 

3.1. coupling in the rotating frame 

Consider the time-dependent linear coupling between two resonator modes of the form 
Hint = — 2Ao sin Udt(a + a)){b + ¥), where Ao is the coupling amplitude and u>d is the 
modulation frequency of the coupling strength. Here the operators a (a') and b (b>) are 
the annihilation (creation) operators for the resonant mode. At the driving frequency 
Ud = uJb — cj a , the coupling in the interaction picture has the form 

Hi = t\ (a r b - aP) (7) 

which generates a beam splitter operation in a similar to fashion to that found in 
quantum optics. At the driving frequency u>a = ou b + u a , the coupling has the form 

Hi = i\ (ab - tilt) (8) 

which generates a squeezing operation and hence entanglement between the two modes. 
These two operations are sufficient for manipulating entanglement and for generating 
squeezed states of the resonators [20J as we will discuss in detail below. 

We focus on the Gaussian states of the coupled resonator modes. A Gaussian state 
can be fully characterized by the covariance matrix: 

with i,j = 1, 2, 3, 4. The operators in this expression are the quadrature variables of the 
resonators with X\ = x a , x 2 = p a , £3 = £&, £4 = Pb for the modes a and b respectively. 
Here () is the ensemble average for the density matrix of the Gaussian states and [, ] is 
the commutator between the variables. Note that linear operations and dissipation due 
to white noise maps Gaussian states to Gaussian states [25J. 

The generic form of the covariance matrix for coupled resonators is 



where A (B, C) is 2 x 2 diagonal matrix with the elements a± and 

A=( a ! ) (111 




(and similarly for 6j and Cj). We set the linear displacement to be (xi) = 0, which 
does not affect the entanglement of the system. To describe the finite temperature of 
the resonators, we use the effective temperature index Q a = coth |^ for a = a,b. 
When k B T huj a , we have Q a = 1 and the mode is in its ground state. The initial 
covariance matrix of the thermal states for the uncoupled resonators is a diagonal matrix 
ao with the diagonal elements (0 Q , a , 0&, 0&)/4 where the quadrature variances of the 
resonators are (xf) = (pf) = 0j/4. 
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3.2. linear operations: squeezing and beam splitter 



It can be shown that the squeezing operation in (JS]) performs the following 
transformation on the covariance matrix, 



a 



where 
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is the squeezing operator and r 2 = A t is the squeezing parameter over a duration t. We 
find that 

/ coshr2 sinhr2 \ 
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Another linear operation, the beam splitter type of operation, on the resonator 
modes in (jHJ), performs the following transformation on the covariance matrix, 
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This transformation produces a beam splitter operation on the coupled resonators. In 
particular, it swaps the states of the two modes when (f> = ir/2. 



4. Entanglement and squeezing of the coupled resonators 



In this section, we study the quantum engineering of entanglement between coupled 
macroscopic resonator modes and the squeezing of the resonators using the parametric 
coupling circuits discussed in the previous sections. 
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Figure 3. (a) The logarithmic negativity En(p) for the two mode squeezed states 
versus the squeezing parameter, (b) The logarithmic negativity Ejy(p) for effective 
polarization entanglement state p p i in (|17p versus displacement parameter «2 with 
ai = 3. In both plots, the upper (blue online) curve is for zero temperature with 
©a = Bb = 1 and the lower (red online) curve is for finite temperature with O a = 4 
and 6b = 1. 



4-1. entanglement 

At zero temperature, the squeezing operation in (jSJ) generates squeezed vacuum states 
between two resonator modes which are also entangled states [26]. Here, we show that 
even at finite temperatures, entanglement can be generated by the squeezing operation 
starting from an initial state with the covariance matrix cr - After a duration t sq with 
the coupling magnitude Ao, the squeezing operation transforms the elements of the 
covariance matrix to 



ci = -c 2 = - (6 a + O b ) sinh \ t sq cosh \ t sq 

which increase with the squeezing parameter r 2 = Xot sq . Note that the above relation 
shows Cj < v 'dibi for % — 1,2 with the equality valid for a large squeezing parameter r 2 . 

To quantitatively characterize the entanglement for the above mixed state, we 
calculate the logarithmic negativity [27] , En(p) = log 2 | \p TA 1 1 , of the partially transposed 
density matrix of this two-mode continuous variable Gaussian state. For Gaussian states, 
the logarithmic negativity can be derived directly from the covariances: 




(14) 
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(15) 



01=1 



where the function F(c a ) is 




for 2c a > 1 

log 2 (2c Q ) for 2c Q < 1 
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and the variables be derived from the roots of the following equation: 

x 4 + 4(det(A) + det(B) - 2 det(C))x 2 + 16 det(cx) = 0. 



The matrices A,B,C were defined previously fTTl for the covariance matrix a [TUl 

In figure [3] (a), we plot E N (p) versus the squeezing parameter r 2 . At zero 
temperature (6 a ,b = 1), the entanglement increases nearly linearly with r 2 . At finite 
temperature, finite squeezing is required to overcome the thermal fluctuations before 
any entanglement can be generated. With large squeezing parameter r 2 > 0.5, the 
logarithmic negativity can be approximated as 



increasing linearly with the squeezing parameter. At r 2 = 0.5, O a = 4 for the 
nanomechanical mode, and 0& = 1 for the electrical mode, we have E^(p) = 0.76 
giving finite entanglement. Hence, at finite temperature, entanglement between coupled 
modes can be generated by applying large squeezing. 

Here, we want to compare the two types of entanglement as apposed to the 
generation of entanglement between resonator modes as described in Ref. [11] using 
a different method. For the case of a resonator coupled to a solid-state qubit, 
parametric pulses can be used to flip the qubit state every half period of the resonator 
mode, n/uj a)b , producing an effective amplification of the resonator displacements [28| [3] . 
For two resonators in a pure state, the entangled state can be expressed as oc 
\ai,a 2 ) + | — a 2 ). This state is comparable to the entangled state between two 
qubits 1 0, 0) + |1, 1) and we call it the effective polarization entangled state. For the 
mixed state at finite temperature, the density matrix of the entangled state can be 
written as 



and we neglect the normalization factor in p v \. The logarithmic negatively for this 
state (figure [3]) has also calculated. After a rapid increase with ct 2 , the logarithmic 
negativity saturates at E^{p) — > 1. At finite temperature, we also have -EW(p) — > 1 



4-2. squeezing 

Following the squeezing operation in f |T41) . entanglement between the resonators can be 
manipulated or adjusted using a subsequent beam splitter operation. Applying the beam 
splitter operation for a duration t — t sq , the covariances can be expressed as, 




(16) 
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Figure 4. The covariances and the logarithmic negativity Em{p) for (a, c) 8 a = 
Ob = 1 and (b, d) O a = 2 and Ob = 1- The dashed line at t = t sq signifies the end 
of the squeezing operation with t sq = 0.5/Ao. Following squeezing, the beam splitter 
operation is applied for tt/Xq. 



C\ = ci cos 2ip H — sin 2<p 

with ip = \ (t — t sq ), and similar relations can be derived for a 2 , b 2 , c 2 . One interesting 
feature is that the covariance matrix is divided into a direct sum of two subsets between 
the variables {x a ,x b } and {p a ,Pb} respectively. Here, the covariances show oscillatory- 
behavior depending on the applied time t — t sq or the phase <p of the applied beam splitter 
operation, as is shown in figure HI 

Starting at maximum entanglement at t = t sq , we see E^{p) decreases to zero near 
t = t sq + 7T/AXQ {ip = 7r/4) when the two resonator modes become separable, then returns 
to its maximum value at t — t sq + 7r/2Ao (<p = 7r/2). At t — t sq + 7r/4Ao, we have 

ai = b 2 = (0 a + 6 )e- 2Ao *-/8 

S 2 = 6i = (e a + fe )e 2Ao '-/8 (19) 

showing a significant amount of squeezing for the quadrature fluctuations a\ = (x^) 
and b 2 = {pi). The cross correlation at this point is c\ = c 2 = (O a — 0&)/8. At 
finite temperature, as plotted in figure H] (d), the negativity decreases to zero and the 
coupled resonators are in a separable state for a finite interval near t = t sq + vr/4A [26J. 
With increasing temperature, the duration of this interval increases as well, where the 
entanglement is diminished by thermal fluctuations. Meanwhile, the correlation between 
the logarithmic negativity and the covariance elements Ci j2 can also be seen in figure H] 
(c) and (d). 

In the special situation of zero temperature with O a = 06, as plotted in figure 0] (a) 
and (c), squeezed states can be generated for the resonator modes. Here we have a\ = b 2 , 
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a 2 = b\ and Ci = — C2 at all times. When t = t sq + 7r/4Ao, the cross correlations as well 
as the entanglement vanish with ci 2 = E^{p) = 0. When both modes are initially 
prepared in their ground states with G flj b = 1, V^i^2 — 1/4. Thus, a squeezed state 
is generated in each resonator. For the coupled system of a nanomechanical resonator 
and an electrical resonator, this provides a novel way of generating squeezed states in 
the nanomechanical mode using linear parametric coupling. 

5. Detection and the covariance matrix 

A crucial requirement for studying the quantum properties of coupled resonators is the 
detection and verification of any engineered entanglement. Below we show that by only 
measuring the quadrature variances of one resonator - the electrical resonator - the full 
covariance matrix can be constructed. 

In reality, measurement of the vibration of a nanomechanical mode is limited by the 
weak coupling between that mode and the detector. The electrical mode, in contrast, 
contains an electromagnetic signal that couples more strongly to the detector than 
the nanomechanical mode. Hence, the electrical mode provides an effective probe of 
the nanomechanical mode, which reduces the demand on the detector efficiency. For 
coupled electrical resonators, both modes can be measured simultaneously to directly 
test the theoretical results presented above. 

The variances of the {x a , Xb} quadratures at any time during the beam splitter 
operation are determined by three initial parameters: a 1; bi, c±, and similarly for the 
{p a ,Pb} quadratures. Measurements of the variances b\ (and 62), at three different times 
(in three sets of experiments with the beam splitter operation applied for different ip) 
that are linearly independent, can provide complete information about the entanglement 
dynamics of the coupled modes. For example, we choose to make measurements at 
ip = 0, 7r/4, tt/2 as defined above. At <p — 0, b\ — b\ is measured; at ip = n/2, b\ = a 1 
is measured; at ip = 7r/4, bi = 01 + 1 + C\ is measured, from which C\ can be obtained 
by combining the previous results. Hence, the variances of the coupled modes can be 
uniquely determined by the measurements of the quadratures of (one of) the electrical 
modes. The times at which the measurements are performed can be adjusted. By 
choosing (p = 7r/8, 7r/4, 37r/8, the quadrature fluctuations exceed thermal noise with 
b\ ^> 0b/4 in all three measurements as is shown in figure H] (a, b). Note that to 
measure the squeezed state of a nanomechanical mode, a fast beam splitter operation 
with a phase it/2 can be applied before measuring the electrical mode. This operation 
transfers the states between the nanomechanical mode and the electrical mode, so that 
the subsequent measurements of the electrical resonator provide direct information of 
the nanomechanical mode. 

This method provides a useful way of detecting the "hard" mechanical mode by 
detecting the "easy" electrical mode. Mechanical modes are in general "hard" to 
measure because they couple very weakly to detectors. By transferring the state of 
the mechanical mode to the electrical mode, which is in general "easy" to detect, we 
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can clearly access the dynamical features of the mechanical mode. This state transfer 
technique was developed previously for phase qubits coupled to an electrical resonator 
[2]. Use of this method yields important information about entanglement by taking 
advantage of a dynamical process to reduce the requirements on the detector efficiency 
For the superconducting transmission line resonator, phase sensitive detection of 
the quadrature variables has been performed to reveal the vacuum Rabi splitting due 
to resonator-qubit coupling by measuring the transmission or reflection of a resonant 
signal [5] and for a resonator coupled to a nanomechanical beam to study the fluctuation 
of the quadrature variables approaching their quantum limit [?,[9]. Alternatively, some 
work has studied the phase sensitive detection of resonator modes with a single electron 
transistor (SET) by mixing the signal of the resonator with a large radio-frequency 
reference signal and coupling the mixed signal to a SET detector [221129] . The nonlinear 
response of the SET enables a high resolution measurement of the electrical resonator 
[30J. The sensitivity of the measurement, however, can be limited by the detector noise. 
In the discussion above, only the detection of the covariance matrix of the coupled 
resonators was studied. Full characterization of the resonator states can be obtained by 
a quantum state tomography method as has been developed in the context of quantum 
optics [3~T] . 

6. Discussion and Conclusions 

The effects studied in this paper can be tested with realistic parameters. Typical 
parameters for a nanomechancial resonator [7j are o; a /27r = 100 MHz, Q a = 10 4 , 
do = 50 nm, and Sxq = 5fm. For a superconducting transmission line resonator 
[5], uj b /2Ti = 5 GHz, cL = 4 fF, and Q b = 10 4 . With a bias voltage of V° = 4 V 
and a coupling capacitance of C° = 0.65 fF, we find Xq/2tt = 6 MHz and the linear 
operations can be performed over a characteristic time scale of t sq = 7r/4Ao ~ 130 ns. For 
two coupled electrical resonators, choosing reasonable parameters for lumped-element 
quantum circuits [32] u a /2n = 9 GHz, u>b/27t = 10 GHz, L a = L b = 5L C = 10M = 
500 pH, I Q = 2.8 /iA, and a modulating flux drive like that described in section [2721 gives 
\q/2ti = 6 MHz so that the linear operations can also be performed for this system over 
a timescale of t sq ~ 130 ns. 

We didn't discuss the effect of decoherence on the dynamics of the coupled resonator 
systems [33]. Many discussions can be found in the literature. The finite quality 
factors of both the electrical resonator and the nanomechanical resonator can limit 
the entanglement or squeezing in this scheme. Here, for a squeezing parameter r 2 = 0.5, 
we have n ~ 2(xl) ~ 1. Assuming a modest quality factor of Q = 10 4 which 
is experimentally realizable, the decoherence time is on the order of 2-nQjifiLOh) ~ 
2 fj,s 3> t sq sufficiently long enough to observe the entanglement and squeezing generated 
in these systems. 

For the measurement process, we proposed a scheme using homodyne detection of 
the covariance matrix of the superconducting electrical resonator mode by applying the 
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beam splitter operation at three different durations giving a full account of the covariance 
matrix of the coupled resonators. This avoids the difficulty of a direct measurement of 
the nanomechancial resonator mode. Considering a coupling capacitance of C° ~ 0.5 fF 
and a maximum possible bias voltage of V x ~ 10 V, the demands on the detector 
efficiency are great if we consider directly measuring the quantum limited fluctuations 
of the nanomechanical resonator. For two coupled electrical resonators, direct detection 
on both resonators together can reveal the covariances for all controlled operations, 
extremely useful for testing these types of systems and the schemes we have described 
for generating entanglement and squeezed states. 

In conclusion, we studied the generation of controllable entanglement and squeezing 
in coupled macroscopic resonators, in particular, a nanomechanical-electrical resonator, 
by applying parametrically modulated linear coupling. Two systems are studied 
in detail: the coupling between a nanomechanical resonator and a superconducting 
transmission line resonator, and the coupling between two superconducting electrical 
resonators. The parametric coupling is calculated in both cases. Furthermore, effective 
detection of the entanglement and squeezing in one of the modes can be achieved 
by measuring the other (electrical) mode. We have considered specific, reasonable 
operating parameters and find that both of these systems should be experimentally 
feasible. Squeezing of the nanomechanical system, although difficult, should be possible 
at very low temperatures whereas the relatively simple two coupled electrical resonator 
system should show clear squeezing for a variety of typical operating conditions. The 
scheme studied here is a key component for entanglement based quantum information 
protocols such as quantum teleportation in a solid-state network [23] . 
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